Lagrangean description of nonlinear dust— ion acoustic waves in dusty plasmas * 
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An analytical model is presented for the description of nonlinear dust-ion-acoustic waves propa- 
gating in an unmagnetized, coUisionless, three component plasma composed of electrons, ions and 
inertial dust grains. The formulation relies on a Lagrangean approach of the plasma fluid model. The 
modulational stability of the wave amplitude is investigated. Difi'erent types of localized envelope 
electrostatic excitations are shown to exist. 
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I. INTRODUCTION 

In the last two decades, dusty plasmas (DP) have at- 
tracted a great deal of attention due to a variety of 
new phenomena observed in them and the novel physical 
mechanisms involved in their description [l|J4|. In addi- 
tion to known plasma electrostatic modes '3|, new oscil- 
latory modes arise in DP 0,0 , amone; which the dust-ion 
acoustic wave (DIAW) and dust acoustic waves (DAW) 
are of significant interest in laboratory dusty plasma dis- 
charges. In the DIAW the restoring force comes from the 
pressures of inertialess electrons, whereas the ion mass 
provides the inertia, similar to the usual ion-acoustic 
waves in an electron- ion plasma. Thus, the DIAW is 
characterized by a phase speed much smaller (larger) 
than the ion (electron) thermal speed, and a frequency 
much higher than the dust plasma frequency uip.d', there- 
fore, on the timescale of our interest, stationary dust 
grains do not participate in the wave dynamics; they just 
affect the equilibrium quasi-neutrality condition. As a 
matter of fact, the DIAW phase velocity is higher than 
that of lA waves, due to the electron density depletion in 
the background plasma when dust grains are negatively 
charged; quite remarkably, this fact results in suppression 
of the Landau damping mechanism 0| , known to prevail 
over the lAW propagation in an electron- ion plasma Q . 

The linear properties of the lAWs have been quite ex- 
tensively studied and now appear well understood 0. 
As far as nonlinear effects are concerned, various stud- 
ies have pointed out the possibility of the formation of 
DIAW-related localized structures, due to a mutual com- 
pensation between nonlinearity and dispersion, including 
small- amplitude pulse solitons, shocks and vortices 
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Furthermore, the propagation of nonlinear ly modulated 
DIA wave packets was studied in Ref. 0, in addition 
to the formation of localized envelope soliton-modulated 
waves due to the modulational instability of the carrier 
waves. A very interesting known approach, not yet in- 
cluded in our current knowledge with respect to the DIA 
plasma waves, is the Lagrangean description of a nonlin- 
ear wave profile. In the context of electrostatic plasma 
waves, this formalism has been employed in studies of 
electron plasma waves 0, 0, Hi an d, more recently, ion- 
acoustic y\ and dust-acoustic "10] waves. Our aim here 
is to extend previous results by applying the Lagrangean 
formalism to the description of nonlinear DIAWs propa- 
gating in dusty plasmas. 

We shall consider the nonlinear propagation of dust- 
ion-acoustic waves in a coUisionless plasma consisting of 
three distinct particle species 'a': an inertial species of 
ions (denoted by 'i'; mass m.^, charge qi — +Zie; e de- 
notes the absolute of the electron charge), surrounded 
by an environment of thermalized electrons (mass rrig, 
charge — e), and massive dust grains (mass M, charge 
Qd = sZdC, both assumed constant for simplicity); Zd de- 
notes the charge state of dust grains; we leave the choice 
of dust grain charge sign s = qd/\qd\ (= — 1/ + 1 for neg- 
ative/positive dust charge) open in the algebra. Charge 
neutrality is assumed at equilibrium. 



II. THE MODEL 

Let us consider the hydrodynamic-Poisson system of 
equations which describe the evolution of the ion 'fluid' 
in the plasma. The ion number density is governed by 
the continuity equation 



— + \/(ni Uj) 



0, 



where the mean velocity obeys 
dvLi „ ZiC 

~dt^ 
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(1) 



(2) 



The electric field E — —V <f> is related to the gradient of 
the wave potential <I>, which is obtained from Poisson's 
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equation \7 ■ 'E = Att qs Ug, viz. 
Alternatively, one may consider 



dE 
'dt 



Qa • 



(3) 



(4) 



We assume a near-Boltzmann distribution for the elec- 
trons, i.e. TT-e « Uefi exp{e^ / ksTf,) (Tg is the electron 
temperature and is Boltzmann's constant). The dust 
distribution is assumed stationary, i.e. rid ~ const.. The 
overall quasi-neutrality condition at equilibrium then 
reads 



Zi riifi + sZd Ud - Ue.o = . 



A. Reduced Eulerian equations 



(5) 



By choosing appropriate physical scales, Eqs. (l)-(3) 
can be cast into a reduced (dimcnsionlcss) form. Let 
us define the ion-acoustic speed Cg — {kBTe/rrii)^^^. 
An appropriate choice for the space and timescales, L 
and T = L/cs, are the effective Debye length A^i = 
{kBTe/AnZfuifie'^y^'^ = Cs/up^i and the ion plasma pe- 
riod ujpl — {AirnifiZfe^ /nii)^^/'^ , respectively. Alter- 
natively, one might leave the choice of L (and thus 
T ~ L/cs) arbitrary - following an idea suggested in 
Refs. 13, ~ which leads to the appearance of a di- 
mensionless dispersion parameter 5 — l/{ujp^iT) = Xd/L 
in the formulae. The specific choice of scale made above 
corresponds to (5 = 1 (implied everywhere in the follow- 
ing, unless otherwise stated); however, we may keep the 
parameter S to 'label' the dispersion term in the forth- 
coming formulae. 

For one-dimensional wave propagation along the x 
axis, Eqs. Q - ® ca.n now be written as 



dn 

W + 
du 

'dt 



d{n u) 

dx 
du 

dx 
dx"^ 



= 0, 



{fi — n) . 



(6) 



where all quantities are dimensionless: n — ni/ni^Q, 
u = Ui/wo and </> = '^/^o] the scaling quantities are, 
respectively: the equilibrium ion density riifi, the ef- 
fective sound speed uo — Cs (defined above) and $0 = 
ksTe/ {Zie). The (reduced) electron and dust back- 
ground density n is defined as 



ZiUio Zini_o 



1 



(7) 



where we have defined the DP parameter fj, = 
ne,o/{Zinifi), and made use of Eq. Note that both 
n and h reduce to unity at equilibrium. 



We shall define, for later reference, the function /(</)) 
n - given by Eq. |(7J) - and its inverse function 



r\x) = Z, In 1 



1 



(8) 



viz. f{(j>) = X implies (f> = f~^{x) = g{x). 

We note that the dependence on the charge sign s is 
now incorporated in /i = 1 + sZdUdfi/ {ZiUi^o); retain 
that /i < 1 (/i > 1) corresponds to negative (positive) 
dust. Remarkably, since the dust-free limit is recovered 
for /i = 1, the results to be obtained in the following are 
also straightforward valid for ion-acoustic waves propa- 
gating in (dust-free) e-i plasma, upon setting /i = 1 in 
the formulae. 

The well-known DIAW dispersion relation uj'^ — 
c^k"^ /{PX^jj -f 1) 3 is obtained from Eqs. (QJ to 
On the other hand, the system yields the reduced re- 
lation uP' = k'^ jip'k^ + 1), which of course immediately 
recovers the former dispersion relation upon restoring di- 
mensions (regardless, in fact, of one's choice of space scale 
L; cf. definition of (5). However, some extra qualitative 
information is admittedly hidden in the latter (dimen- 
sionless) relation. Should one consider a very long space 
scale L ^ X^ (i.e. (5 <C 1), one readily obtains uj ^ k 
(unveiling the role of (5 as a characteristic dispersion con- 
trol parameter). Finally, the opposite limit of short L 
(or infinite S) corresponds to ion plasma oscillations (viz. 
Lu = ujp,i ~ constant). 



B. Lagrangean description 

Let us introduce the Lagrangean variables 
which are related to the Eulerian ones {x, t\ via 



S, = X - I u{i, t') dr , 
Jo 



(9) 



See that they coincide at i = 0. Accordingly, the space 
and time gradients are transformed as 

d/dx -> d/d^ , d/dt ^ d/dr - u d/d^ , 

where we have defined the quantity 



(10) 



Note that the convective derivative D = d/dt + ud/dx 
is now plainly identified to d/dr. Also notice that a 
satisfies a(^, t = 0) — and 



da{C,T) du{^,T) 



dr 



di 



(11) 



As a matter of fact, the Lagrangean transformation de- 
fined here reduces to a Galilean transformation if one 
suppresses the evolution of w, i.e. for u = const, (or 
du/dr = du/d^ = 0, hence a = 1). Furthermore, if 
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one also suppresses the dependence in time r, this trans- 
formation is reminiscent of the travelling wave ansatz 
f{x,t) = f{x — vt = s), which is widely used in the 
Sagdeev potential formalism [3. ITU. 

The Lagrangean variable transformation defined above 
leads to a new set of reduced equations 



a 



dr 



(e,T)n(C,0) 



(12) 



a 

-—r (^'-)^^(13) 

47rZie[n(^,r) - nrii^o] (14) 
-47rZ,en(e, r)M(^, r) , (15) 



where we have temporarily restored dimensions for physi- 
cal transparency; recall that the (dimensionless) quantity 
n, which is in fact a function of 0, is given by Q. One 
immediately recognizes the role of the (inverse of the) 
function «(£, r) as a density time evolution operator; cf. 
Eq. \i2\ (12|| . Poisson's equation is now obtained by 
eliminating (p from Eqs. (|13l I14|l 



nn^^o) 



(16) 



Note that a factor 5"^ should appear in the left-hand side 
if one rescaled Eq. as described above; cf. the last 
of Eqs. This will be retained for later reference, 

with respect to the treatment suggested in Ref. 9] (see 
discussion below). 

In principle, our aim is to solve the system of Eqs. 
(O to lO or, by ehminating cf), Eqs. (jHl), (O and 
(|16|l for a given initial condition n(^,r = 0) = no(0, 
and then make use of the definition ^ in order to invert 
back to the Eulerian arguments of the state moment vari- 
ables (i.e. density, velocity etc.). However, this abstract 
scheme is definitely not a trivial task to accomplish. 



III. NONLINEAR DUST-ION ACOUSTIC 
OSCILLATIONS 

Multiplying Eq. by m(5, t) and then adding to Eq. 
(|15|l . one obtains 



dE{£„T) 
dr 



Combining with Eq. 113|l . one obtains 



(17) 



(18) 



where LOp^i is the ion plasma frequency (defined above). 
Despite its apparent simplicity, Eq. H18|) is neither an 



ordinary differential equation (ODE) - since all variables 
depend on both time r and space ^ - nor a closed evo- 
lution equation for the mean velocity u(^,t): note that 
the (normalized) background particle density n depends 
on the potential (/) and on the plasma parameters; see its 
definition (0. The evolution of the potential 0(^,t), in 
turn, involves u{^,t) (via the quantity a(f,r)) and the 
ion density n(^,T). 

Eq. H18(l suggests that the system performs nonlinear 
oscillations at a frequency uj — Upiv}/"^. Near equilib- 
rium, the quantity fi is approximately equal to unity and 
one plainly recovers a linear oscillation at the ion plasma 
frequency ujp^i- Quite unfortunately this apparent sim- 
plicity, which might in principle enable one to solve for 
m(^, r) and then obtain {^, r} in terms of {x, i} and vice 
versa (cf. Davidson's treatment for electron plasma os- 
cillations in Ref. 0; also compare to Ref. Q], setting 
7 = therein), is absent in the general (off-equilibrium) 
case where the plasma oscillations described by Eq. H18|l 
are intrinsically nonlinear. 

Since Eq. H18|l is in general not a closed equation for 
M, unless the background density n is constant (i.e. in- 
dependent of <j), as in Refs. 7, 8.]), one can neither apply 
standard methods involved in the description of nonlin- 
ear oscillators on Eq. ifT^ (cf. Ref. 0), nor reduce the 
description to a study of Eqs. 1)181 I17|l (cf. Ref. '^), 
but rather has to retain all (or rather five) of the evolu- 
tion equations derived above, since five inter-dependcnt 
dynamical state variables (i.e. n, u, i?, <j) and a) are in- 
volved. This procedure will be exposed in the following 
Section. 



IV. PERTURBATIVE NONLINEAR 
LAGRANGEAN TREATMENT 

Let us consider weakly nonlinear oscillations performed 
by our system close to (but not at) equilibrium. The basis 
of our study will be the reduced system of equations 



d_ 

dr 



[an] 



0, 



du 

OT 

dE 

— = {n-n)a, 

da du 
9^ " 9C ' 



(19) 



which follow from the Lagrangean Eqs. (|12() to l|16|l by 
scaling over appropriate quantities, as described in ^11 Al 
[l3l |. This system describes the evolution of the state 
vector, say S = {a,n,u,E,(j)) (g 3?^), in the Lagrangean 
coordinates defined above. We will consider small devi- 
ations from the equilibrium state So = (1, 1, 0, 0, 0)"^, by 



taking S = S(0) -I- eSi 



(0) 



,(0) 



where e (^ 1) 
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is a smallness parameter. Accordingly, we shall Tay- 
lor develop the quantity n{(j)) near « 0, viz. (/) « 

,^ , ,2 J. , • ^ _„J_„ i_ ~ 



£01 + 



in order to express n as 

/,2 , „ -3 



- ci0 4- 
eci0i 



+ C30 



C20?) 



3(C103 +2C20102 +C30?) 



(20) 



where the coefficients Cj (j = 1,2,...), which are deter- 
mined from the definition of h, contain all the es- 
sential dependence on the plasma parameters, e.g. /x; 
making use of ~ X^i^Lo x"^ /i^^-i one readily obtains 



ci = ^i/Zi 



C2 = /i/(2Z2) , C2 - ^l/i6Zf) . 



Remember that for = 1 (i.e. for vanishing dust) one 
recovers the expressions for lAWs in e-i plasma. 

Following the standard reductive perturbation tech- 
nique we shall consider the stretched (slow) La- 
grangean coordinates Z = e(^— Ar), T = t (where 
A G 3? will be determined later). The perturbed state of 
(the j— th — J = 1, 5 — component of) the state vec- 
tor S'"' is assumed to depend on the fast scales via the 
carrier phase 9 = k^ — cot, while the slow scales enter the 
argument of the {j—th element's) th harmonic ampli- 



tude S. 



3-1 



(where sj"^; 



'• Sj 



(")* 



EZ-o. S^{Z, T)e 



_( -iji ensures reality), 
derivative operators as 



Treating the 



d_ 



d_ 



eA 



d_ 

dZ 



d_ 

df 



d_ 



d_ 



dZ ' 



and substituting into the system of evolution equations, 
one obtains an infinite series in both (perturbation or- 
der) e" and (phase harmonic) I. The standard pertur- 
bation procedure now consists in solving in successive 
orders ~ e" and substituting in subsequent orders. The 
method involves a tedious calculation which is however 
straightforward; the details of the method are presented 
e.g. in Ref. so only the essential stepstones need to 
be provided here. 

The equations obtained for n ~ I — 1 determine the 
first harmonics of the perturbation 



.(1) - 



E 



(fcV^')^, 



(1) 



(21) 



where ij) denotes the potential correction (f)\ . The 
cyclic frequency to obeys the dispersion relation uj"^ — 
k'^ / {k'^ + sci) , which exactly recovers, once dimensions 
are restored, the standard lAW dispersion relation ^ 
mentioned above. 

Proceeding in the same manner, we obtain the second 
order quantities, namely the amplitudes of the second 

(2) (2) 

harmonics S2 and constant ('direct current') terms Sg , 

(2) 

as well as a finite contribution Sj^ to the first harmonics; 
as expected from similar studies, these three (sets of 5, at 



each n, /) quantities are found to be proportional to ^/i^, 
IV'P and d-ip/dZ respectively; the lengthy expressions are 
omitted here for brevity. The (71 = 2, Z = 1) equations 
provide the compatibility condition: A = lo{1 — u}'^)/k = 
duj / dk; A is therefore the group velocity Vg{k) = uj'{k) at 
which the wave envelope propagates. It turns out that Vg 
decreases with increasing wave number fc; nevertheless, 
it always remains positive. 

In order ~ e'^, the equations for I = 1 yield an ex- 
plicit compatibility condition in the form of a nonlinear 
Schrodinger-type equation (NLSE) 



* dT 

^(1) 



+ g|V^|2^/. = o. 



(22) 



Recall that tp = (f>i denotes the amplitude of the first- 
order electric potential perturbation. The 'slow' variables 
{Z,T} were defined above. 

The dispersion coefficient P is related to the curvature 
of the dispersion curve as P = uj"{k)/2 = — 3ti;^(l — 
Lu'^) / {2k'^) . One may easily check that P is negative (for 
all values of fc). 

The nonlinearity coefficient Q is due to carrier wave 
self-interaction. It is given by the expression 



12 fc4 Zf 



3Zfk' 



3 7.6 



'3{n+4)Zfk'^+3(l-2fi-5n^)Z,k'^-fi{3fi- 

(23) 



where the coefficients ci.2,3 were defined above. 
For low wave number fc, Q goes to —00 as 



i)V/' 1 



5/2 



12 z. 



A. Modulational stability analysis 

According to the standard analysis we can lin- 
earize around the plane wave solution of the NLSE H22|l 
■0 — ipe^QW -(- c.c. , (c.c: complex conjugate) - no- 
tice the amplitude dependence of the frequency shift 
Aw = e^QIV'P - by setting tp = ^po + eipi , and then 
assuming the perturbation ipi to be of the form: ijji = 
c.c. Substituting into H22() . one thus 



^^^^ ^rikC-ur) 

readily obtains a)^ 



P2fc2 p _ 2(Q/P)|0i^or • The 



wave will thus be stable (V k) if the product PQ is nega- 
tive. However, for positive PQ > 0, instability sets in for 
wavenumbers below a critical value k^r — y^2Q/P |'0i,o|j 
i.e. for wavelengths above a threshold Ac- — 2'K/kcr] 
defining the instability growth rate cr = \I'ml>{k)\, we see 
that it reaches its maximum value for k = fccr/\/2, viz. 

(y-max = \lTT^'^\k=k,,/./2 ^ 1*51 IV'l.ol^ • 

We see that the instability condition depends only on the 
sign of the product PQ, which may be studied numeri- 
cally, relying on the expressions derived above. 
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B. Finite amplitude nonlinear excitations 

The NLSE ^ is long known to possess distinct types 
of localized constant profile (solitary wave) solutions, de- 
pending on the sign of the product PQ |ElillI3- Re- 
member that this equation here describes the evolution 
of the wave's envelope, so these solutions represent slowly 
varying localized envelope structures, confining the (fast) 
carrier wave. The analytic form of these excitation can be 
found in the literature (see e.g. in for a brief review) 
and need not be derived here in detail. Let us however 
briefly summarize those results. 

Following Ref. , we may seek a solution of Eq. H22|l 
in the form 7/;(C, r) = p{Z, T) e'®'-^''^''-|-c.c., where p, a are 
real variables which are determined by substituting into 
the NLSE and separating real and imaginary parts. The 
different types of solution thus obtained arc summarized 
in the following. 

For PQ > we find the (bright) envelope soliton 

p^±posech(^^^^y Q^l^[u,Z-{n+^ul)T], 

(24) 

which represents a localized pulse travelling at the enve- 
lope speed Ue and oscillating at a frequency Q (at rest). 
The pulse width L depends on the maximum amplitude 
square po as L = (2P/(5)^^^/po- Since the product PQ 
is always positive for long wavelengths, as we saw above, 
this type of excitation will be rather privileged in dusty 
plasmas. The bright-type envelope soliton is depicted in 

Fig. nil, b. 

For PQ < 0, we obtain the dark envelope soliton (hole) 

M 



1 — sech 



-,1/2 



±pi tanh 



Z — UeT 

L' 



e = — 



1 

2P 



u„Z -\\ul- 2PQp^ ) T 



(25) 



which represents a localized region of negative wave den- 
sity (shock) travelling at a speed Ue] see Fig. Again, 
the pulse width depends on the maximum amplitude 
square pi via L' = (2|P/Q|)1/Vpi- 

Finally, still for PQ < 0, one also obtains the gray 
envelope solitary wave 



1 — a sech 



Z — u^T 



1/2 



(26) 



which also represents a localized region of negative wave 
density. Comparing to the dark soliton H25|) . we note 
that the maximum amplitude p2 is now finite (non-zero) 
everywhere; see Fig. The the pulse width of this 

gray-type excitation L" — yj2\P/Q\/{a P2) now also de- 
pends on an independent parameter a which represents 



the modulation depth (0 < a < 1). The lengthy expres- 
sions which determine the phase shift Q and the param- 
eter a, which are omitted here for brevity, can be found 
in Refs. For a = 1, one recovers the dark soliton 

presented above. 



An important qualitative result to be retained is that 
the envelope soliton width L and maximum amplitude 
p satisfy Lp ^ \fPlQ (see above), and thus depend on 
(the ratio of) the coefficients P and Q\ for instance, re- 
gions with higher values of P (or lower values of Q) will 
support wider (spatially more extended) localized exci- 
tations, for a given value of the maximum amplitude. 
Contrary to the KdV soliton picture, the width of these 
excitations does not depend on their velocity. It does, 
however, depend on the plasma parameters, e.g. here p. 



The localized envelope excitations presented above 
represent the slowly varying envelope which confines the 
(fast) carrier space and time oscillations, viz. — 
^'(X, Z) cos(fc^ — lot) for the electric potential (and 
analogous expressions for the density rii etc.; cf. (j2B). 
The qualitative characteristics (width, amplitude) of 
these excitations, may be investigated by a numerical 
study of the ratio P/Q = ri{k; p): recall that its sign de- 
termines the type (bright or dark) of the excitation, while 
its (absolute) value determines its width for a given am- 
plitude (and vice versa). In Fig. |21we have depicted the 
behaviour of as a function of the wavenumber k and 
the parameter p: higher values of p correspond to lower 
curves. Remember that, for any given wavenumber k, 
the dust concentration (expressed via the value of p) de- 
termines the soliton width L (for a given amplitude p; 
see discussion above) since L ~ Therefore, we 

see that the addition of negative dust generally (p < 1) 
results to higher values of 77 (i.e. wider or higher soli- 
tons), while positive dust (/i > 1) has the opposite effect: 
it reduces the value of rj (leading to narrower or shorter 
solitons). In a rather general manner, bright type soli- 
tons (pulses) seem to be rather privileged, since the ratio 
■q (or the product PQ) of the coefficients P and Q is pos- 
itive in most of the k,p plane of values. One exception 
seems to be very the region of very low values of p (typ- 
ically below 0.2), which develops a negative tail of r/ for 
small k (< 0.3 kn): thus, a very high (> 80 per cent) elec- 
tron depletion results in pulse destabilization in favour of 
dark-type excitations (Fig. d). Strictly speaking, 77 
also becomes negative for very high wave number values 
(> 2.5 fcu); nevertheless, we neglect - for rigor - this re- 
gion from the analysis, in this (long wavelength A) fluid 
picture (for a weak dust presence, short A DIAWs may 
be quite strongly damped; however, this result may still 
be interesting for a strong presence of dust, when Landau 
damping is not a significant issue 1]). 
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V. RELATION TO PREVIOUS WORKS: AN 
APPROXIMATE NONLINEAR LAGRANGEAN 
TREATMENT 

By combining the Lagrangean system of Eqs. H12() to 
one obtains the (reduced) evolution equation 



no 



n d 
no 



where the function g(x) [defined in Eq. 
at 



(27) 

is evaluated 



w{n) = n 



Note that the ion density n has been scaled by its equi- 
librium value ni o, to be distinguished from the initial 
condition no = n(^, t — 0). 

Despite its complex form, the nonlinear evolution 
equation (|?7jl can be solved exactly by considering dif- 
ferent special cases, as regards the order of magnitude 
of the dispersion-related parameter S. This treatment, 
based on Ref. Q, will only be briefly summarized here, 
for the sake of reference. 

First, one may consider very short scale variations, i.e. 
L X^i {or S ::§> 1). This amounts to neglecting collec- 
tive effects, so oscillatory motion within a Debye sphere is 
essentially decoupled from neighboring ones. By consid- 
ering w{n) w — (5^ fi d'^ (l/n) / dr'^ and (p ^ (i.e. n w 1), 
one may combine Eqs. (|15^ and H27|) into 



9r2 



--11=0. 

n 



(28) 



which, imposing the initial condition n{^,0) = no(^), 
yields the solution 



"(e,r) 



noiO 



lioii) +(l_!ioii))cost^p,T 



(29) 



Note that if the system is initially at equilibrium, viz. 
no(^) = nifl, then it remains so at all times r > 0. Now, 
one may go back to Eq. ifT^ and solve for a(f , r), which 
in turn immediately provides the mean fluid velocity u 

u(^,t) = Ljp iSinujpsT [ (l - ^ 

via Hl()(l . and then E{^,t) and (/)(^,t). Finally, the vari- 
able transformation (^J may now be inverted, immedi- 
ately providing the Eulerian position x in terms of ^ and 
T. We shall not go into further details regarding this pro- 
cedure, which is essentially analogue (yet not identical) 
to Davidson's treatment of electron plasma oscillations. 

Quite interestingly, upon neglecting the dispersive ef- 
fects, i.e. setting 5 = 0, Eq. ifTTI) may be solved by 
separation of variables, and thus shown to possess a 
nonlinear special solution in the form of a product, say 



iT-i^iT) = ni{0n2{T) This calculation was put for- 

ward in Ref. (where the study of lAW - in a single 
electron temperature plasma - was argued to rely on an 
equation quasi- identical to Eq. (|27|) '). However, the so- 
lution thus obtained relies on doubtful physical grounds, 
since the assumption (5 ~ 0, which amounts to remaining 
close to equilibrium - cf. the last of Eqs. ©, implies an 
infinite space scale L (recall the definition of S) , contrary 
to the very nature of the (localized) nonlinear excitation 
itself. Rather not surprisingly, this solution was shown 
in Ref. Q to decay fast in time, in both Eulerian and 
Lagrangean coordinates. Therefore, we shall not pursue 
this analysis any further. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied the nonlinear propagation of dust ion 
acoustic waves propagating in a dusty plasma. By em- 
ploying a Lagrangean formalism, we have investigated 
the modulational stability of the amplitude of the propa- 
gating dust ion acoustic oscillations and have shown that 
these electrostatic waves may become unstable, due to 
self interaction of the carrier wave. This instability may 
either lead to wave collapse or to wave energy localiza- 
tion, in the form of propagating localized envelope struc- 
tures. We have provided an exact set of analytical ex- 
pressions for these localized excitations. 

This study complements similar investigations which 
relied on an Eulerian formulation of the dusty plasma 
fluid model ^5J. In fact, the Lagrangean picture pro- 
vides a strongly modified nonlinear stability profile for 
the wave amplitude, with respect to the previous (Eule- 
rian) description; this was intuitively expected, since the 
passing to Lagrangean variables involves an inherently 
nonlinear transformation, which inevitably modifies the 
nonlinear evolution profile of the system described. How- 
ever, the general qualitative result remains in tact: the 
dust ion acoustic-type electrostatic plasma waves may 
propagate in the form of localized envelope excitations, 
which are formed as a result of the mutual balance be- 
tween dispersion and nonlinearity in the plasma fluid. 
More sophisticated descriptions, incorporating e.g. ther- 
mal or collisional effects, may be elaborated in order to 
refine the parameter range of the problem, and may be 
reported later. 
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Figure Captions 



Figure 1. 

A heuristic representation of wave packets modulated 
by solutions of the NLS equation. These envelope excita- 
tions are of the: (a, b) bright type {PQ > 0, pulses); (c) 
dark type, (d) gray type {PQ < 0, voids). Notice that 
the amplitude never reaches zero in (d). 

Figure 2. 

The ratio rj = P/Q oi the coefficients in the NLSE 
(|22|l is depicted versus the wave number k (normalized 
over ki)), for several values of the dust parameter fi; in 
descending order (from top to bottom): 0.8, 0.9, 1.0, 1.1, 
1.2. 
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